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The readout of a classical memory can be modelled as a problem of quantum channel discrimination, where 
a decoder retrieves information by distinguishing the different quantum channels encoded in each cell of the 
memory [S. Pirandola, Phys. Rev. Lett. 106, 090504 (2011)]. In the case of optical memories, such as 
CDs and DVDs, this discrimination involves lossy bosonic channels and can be remarkably boosted by the 
use of nonclassical light (quantum reading). Here we generalize these concepts by extending the model of 
memory from single-cell to multi-cell encoding. In general, information is stored in a block of cells by using 
a channel-codeword, i.e., a sequence of channels chosen according to a classical code. Correspondingly, the 
readout of data is realized by a process of "parallel" channel discrimination, where the entire block of cells is 
probed simultaneously and decoded via an optimal collective measurement. In the limit of an infinite block 
we define the quantum reading capacity of the memory, quantifying the maximum number of readable bits per 
cell. This notion of capacity is nontrivial when we suitably constrain the physical resources of the decoder. 
For optical memories (encoding bosonic channels), such a constraint is energetic and corresponds to fixing the 
mean total number of photons per cell. In this case, we are able to prove a separation between the quantum 
reading capacity and the maximum information rate achievable by classical transmitters, i.e., arbitrary classical 
mixtures of coherent states. In fact, we can easily construct nonclassical transmitters that are able to outperform 
any classical transmitter, thus showing that the advantages of quantum reading persist in the optimal multi-cell 
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I. INTRODUCTION 



One of the central problems in the field of quantum infor- 
mation is the statistical discrimination of quantum states [fil- 
ial]. This is a fundamental issue in many protocols, including 
those of quantum communication ||5rlal and quantum cryp- 
tography I9I 4I3II . A similar problem is the statistical discrim- 
ination of quantum channels, also called "quantum channel 
discrimination" (QCD) 11411 . In its basic formulation, QCD 
involves a discrete ensemble of quantum channels which are 
associated with some a priori probabilities. A channel is ran- 
domly extracted from the ensemble and given to a party who 
tries to identify it by using input states and output measure- 
ments. The optimal performance is quantified by a mini- 
mum error probability which is generally non-zero in the pres- 
ence of constraints (e.g., for fixed number of queries or re- 
stricted space of the input states). In general, this is a double- 
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optimization problem whose optimal choices are unknown, 
a feature which makes its exploration non-trivial. Moreover 
QCD may also involve continuous ensembles. A special case 
is the "quantum channel estimation" where the ensemble is in- 
dexed by a continuous parameter with flat distribution. Here 
the goal is to estimate the unknown parameter with minimal 
uncertainty ifrHflrjll . 

Besides its difficult theoretical resolution, QCD is also in- 
teresting for its potential practical implementations. For in- 
stance, it is at the basis of the decoding procedure of the two- 
way quantum cryptography IU7I1 where the secret information 
is encoded in a Gaussian ensemble of phase-space displace- 
ments. Furthermore QCD appears also in the quantum illumi- 
nation of targets IU8l[l9ll . where the sensing of a remote low- 
reflectivity object in a bright thermal environment corresponds 
to the binary discrimination between a very noisy/lossy chan- 
nel (presence of target) and a completely depolarizing channel 
(absence of target). 

More recently, QCD has been connected with another fun- 
damental task: the readout of classical digital memories 112011 - 
Thanks to this connection, Ref. 12011 has laid the basic ideas 
of treating digital memories, such as optical disks, in the field 



of quantum information theory (see also the following stud- 
ies of Refs. lulMZl ). The storage of data, i.e., the writing 
of the memory, corresponds to a process of channel encod- 
ing, where information is recorded into a cell by storing a 
quantum channel picked from some pre-established ensem- 
ble. Then the process of readout corresponds to the process 
of channel decoding, which is equivalent to discriminate be- 
tween the various channels of the ensemble. This is done by 
probing the cell using an input state, also called "transmitter", 
and measuring the output by a suitable detector or "receiver"'. 
Ref. 02OI1 developed this model directly in the bosonic setting, 
in order to apply the results to optical memories, such as CDs 
and DVDs. The central investigation regarded the comparison 
between classical and nonclassical transmitters, where "clas- 
sical transmitters" correspond to probabilistic mixtures of co- 
herent states and encompass all the sources of light which are 
used in today's data storage technology. By contrast, "non- 
classical transmitters" are only produced in quantum optics' 
labs and the y ar e typically based on entangled, squeezed or 
Fock states J24l 12511 . As shown by Ref. ll20Tl . we can con- 
tract nonclassical transmitters that are able to outperform any 
classical transmitter. In particular, this happens in the regime 
of low energy, where a few photons are irradiated over each 
cell of the memory. This regime is particularly important for 
the non trivial implications it can have in terms of increas- 
ing data-transfer rates and storage capacities. Following the 
terminology of Ref. luOll . we call "quantum reading" the use 
of nonclassical transmitters to read data from classical digital 
memories. 

The main results on the quantum reading of memories re- 
garded the single-cell scenario, where each memory cell is 
written and read independently from the others. However, a 
supplementary analysis of Ref. 12011 also showed that the ad- 
vantages of quantum reading persist when we extend the en- 
coding of information from a single- to a multi-cell model. 
Assuming a block-encoding of data, one can use error cor- 
recting codes which make the readout flawless. In this sce- 
nario it is possible to show that the error correction overhead 
can be made negligible at low energies only when we adopt 
nonclassical transmitters B26I1 . Motivated by this analysis, the 
present work provides a full general treatment of the quantum 
reading of memories in the multi-cell scenario. This is done 
by formalizing the most general kind of classical digital mem- 
ory. In this model, information is stored in a block of cells by 
using a channel-codeword, i.e., a sequence of channels cho- 
sen according to some classical code. Then, the readout of 
data is realized by a process of parallel channel discrimina- 
tion. This means that the entire block of cells is probed in 
parallel and then decoded by an optimal collective measure- 
ment. Such a description encompasses all the possible en- 
coding and decoding strategies. Since the storage capacity of 
classical memories is usually very large, an average memory 
is made by a large number of these encoding blocks. The opti- 
mal scenario corresponds to the case where the whole memory 
is represented by a single, very large, encoding block which 
is read in a parallel fashion. In this limit (infinite block) we 
can provide a simple characterization of the memory and re- 
sort to the Holevo bound to quantify the amount of readable 



information. This enables us to define the quantum reading 
capacity of the classical memory, which corresponds to the 
maximum readable information per cell. If we do not impose 
constraints, this capacity equals exactly the amount of infor- 
mation stored in each cell of the memory. However, this is no 
longer the case when we introduce physical constraints on the 
resources accessible to the reading device. In the case of op- 
tical memories, which involve the discrimination of bosonic 
channels, the energy constraint is the most fundamental [5]. 
Thus the quantum reading capacity is properly formulated for 
fixed input energy. This means that we fix the mean total num- 
ber of photons irradiated over each cell of the memory. The 
computation of this capacity would be very important at the 
low energy regime, which is the most interesting for its poten- 
tial implications. Despite its calculation is extremely difficult, 
we are able to provide lower bounds for the most basic op- 
tical memories, i.e., the ones based on the binary encoding 
of lossy channels. For these memories we are able to derive 
a simple lower bound which quantifies the maximum infor- 
mation readable by classical transmitters. We call this bound 
the "classical reading capacity" of the memory and represents 
an extension to the multi-cell scenario of the "classical dis- 
crimination bound" introduced in Ref. 12011 . Remarkably, the 
optimal classical transmitter which irradiates n mean photons 
per cell can be realized by using a single coherent state with 
the same mean number of photons. Thanks to this result, we 
can easily investigate if a particular nonclassical transmitter 
is able to outperform any classical transmitter. This is indeed 
what we find in the regime of few photons. Thus, in the low 
energy regime, we can prove the separation between the quan- 
tum reading capacity and the classical reading capacity, which 
is equivalent to state that the advantages of quantum reading 
persist in the optimal multicell scenario. 



The paper is organized as follows. In Sees. [Hi and Hill we 
review some of the key -points of Ref. J20I1 and its supplemen- 
tary materials, which are preliminary for the new results of 
Sees. IIVWIII In particular, in Sec. HU we review the basic 
notions regarding the memory model with single-cell encod- 
ing. Then, in the following Sec. [Ill] we discuss the simplest 
example of optical memory and its quantum reading. Once 
we have reviewed these notions, we introduce the model with 
multi-cell encoding in Sec. [IV] In Sec. [V] we take the limit 
for infinite block size and we define the quantum reading ca- 
pacity of the memory, both unconstrained and constrained. In 
particular, we specialize the constrained capacity to the case 
of optical memories (bosonic channels). In Sec. [VI] we com- 
pute the lower bound relative to classical transmitters, i.e., the 
classical reading capacity. In the following Sec. I VIII we prove 
that this bound is separated, by showing simple examples of 
nonclassical transmitters which outperform classical ones in 
the regime of few photons. Finally, Sec. I Villi is for conclu- 
sions. 



II. BASIC MODEL OF MEMORY: SINGLE-CELL 
ENCODING 

In the more abstract sense, a classical digital memory can 
be modelled as a one-dimensional array of cells (the general- 
ization to two or more dimensions is just a matter of technical- 
ities). The writing of information by some device or encoder, 
that we just call "Alice" for simplicity, can be modelled as 
a process of channel encoding 12011 . This means that Alice 
has a classical random variable X = {x,p x } with k values 
x = 0, ■ ■ • , k — 1 distributed according to a probability dis- 
tribution p x . Each value x is then associated with a quantum 
channel cf> x via one-to-one correspondence 



X -H> 



thus defining an ensemble of quantum channels 

$ = {<t> x ,p x } • 



(1) 



(2) 



Mathematically speaking, each channel of the ensemble is a 
completely positive trace-preserving (CPT) map acting on the 
state space T>(H) of some chosen quantum system (Hilbert 
space T-L). Furthermore, the various channels are different 
from each other. This means that, for any pair (j> x and (f> x i, 
there is at least one state p <G T>(H) such that 



F[cj> x {p),4> x ,{p)}<l 



(3) 



where Ftp, a) — \^{^/pcr^/p) 1 ^ 2 ] 2 is the quantum fi- 
delity 12711 . Thus, in order to write information, Alice ran- 
domly picks a quantum channel <\> x from the ensemble and 
stores it in a target cell. This operation is repeated identically 
and independently for all the cells of the memory, so that we 
can characterize both the cell and the memory by specifying 
$ (see Fig. [TJ. 

The readout of information corresponds to the inverse pro- 
cess, which is channel decoding or discrimination. The writ- 
ten memory is passed to a decoder, that we call "Bob", who 
queries the cells of the memory one by one. To retrieve in- 
formation from a target cell, Bob exploits a transmitter and 
a receiver. In the simplest case this means that Bob inputs a 
suitable quantum state p and measures the corresponding out- 
put state p x = <p x (p) recording the specific quantum channel 
stored in that cell (see Fig. [T). Note that, given some input 
state p, the ensemble of the possible output states {<fi x (p) , p x } 
is generally made by non-orthogonal states which, therefore, 
cannot be perfectly distinguished by a quantum measurement. 
In other words, the discrimination cannot be perfect and the 
quantum detection will output the correct value x up to an er- 
ror probability P err . It is clear that the main goal for Bob 
is to optimize input state and output measurement in order to 
retrieve the maximal information from the cell. 



A. Multi-copy probing and optical memories 

In a classical digital memory information is stored 
(quasi)permanently. This means that the association between 
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FIG. 1: Basic process of storage and readout. A memory cell 
can be characterized by an ensemble of quantum channels $ = 
{<fix,Px}- Alice picks a quantum channel <f> x (with probability p x ) 
and stores it in a target cell. In order to read the information, Bob 
exploits a transmitter and a receiver. In the simplest scenario, this 
corresponds to inputting a suitable quantum state p and measuring 
the output p x = cj)x(p) by a suitable detector. The detector gives 
the correct answer x up to some error probability P err . Multi-copy 
probing. Since the cell encodes the quantum channel in a stable way, 
we can probe the cell many times. This means that, more gener- 
ally, Bob can input a multipartite state p(s) 6 T>(H® S ) which de- 
scribes s quantum systems. As a consequence, the output will be 
px(s) = </>f s [p(s)], whose global detection gives x up to an error 
probability (which is non-increasing in s). Optical memory. The en- 
coded channel <f> x is a bosonic channel (in particular, single-mode). 
In this case, Bob uses an input state p(s, n) describing s bosonic 
modes and irradiating n mean photons over the cell. 



a single cell and the channel-encoding $ must be stable. As 
a result, Bob can probe the cell many times by using an in- 
put state living in a bigger state space. Given some quantum 
channel 



<f> x : V(H) -»■ V(H) 



(4) 



Bob can input a multipartite state p(s) <G 2?("H® S ) with integer 
a > 1, i.e., describing s quantum systems. As a result, the 
output state will be 



Px ( S ) = ^ripis)} 



(5) 



This state is detected by a quantum measurement applied to 
the whole set of s quantum systems (see Fig. [T). Physically, 
if we consider the process in the time domain, p(s) describes 
the global state of s systems which are sequentially transmit- 
ted through the cell. In other words, the number s can also 
be regarded as a dimensionless readout time [20]. Intuitively, 
it is expected that the optimal P err is a decaying function of 
s, so that it is always possible to retrieve all the information 
in the limit for s — > oo. This suggests that the readout prob- 
lem is nontrivial only if we impose constraints on the physical 
resources that are used to probe the memory. In the case of 
discrete variables (i.e., finite-dimensional Hilbert space) the 
constraint can be stated in terms of fixed or maximum readout 
time s. 

More fundamental constraints come into play when we 
consider an optical memory, which can be defined as clas- 
sical memory encoding an ensemble $ of bosonic channels. 



In particular, these channels can be assumed to be single- 
mode. Since the underlying Hilbert space is infinite in the 
bosonic setting, one has unbounded operators such as the en- 
ergy. Clearly, if we allow the energy to go to infinite, the 
discrimination of (different) bosonic channels is always per- 
fect. As a result, the readout of optical memories has to be 
modelled as a channel discrimination problem where we fix 
the input energy. The minimal energy constraint corresponds 
to fixing the mean total number of photons n irradiated over 
each memory cell 112011 . Thus, for fixed n, the aim of Bob is 
to optimize input (i.e., number of bosonic systems s and their 
state p) and the output measurement. In the following we ex- 
plicitly formalize this constrained problem. 

Let us consider an optical memory with cell $ = {<fi x ,p x } 
where each element (\> x is a single-mode bosonic channel. 
Then, we denote by p(s, n) a multimode bosonic state p e 
T>(H® S ) with mean total energy Ti(pn) = n, where n is the 
total number operator over T-L® s . In other words, this state de- 
scribes s bosonic systems which irradiate a total of n mean 
photons over the target cell (see also Fig. [TJ. We refer to 
the pair (s, n) as to the signal profile. In the bosonic setting 
the parameter s can be interpreted not only as the number of 
temporal modes (therefore, readout time) but equivalently as 
the number of frequency modes, thus quantifying the "band- 
width" of the signal 12011 , Now, for a given input p = p(s, n) 
to the cell $, we have the output state 



p x (s,n) = <t>rip(s,n)} 



(6) 



This output is subject to a quantum measurement over the s 
modes which is generally described by a positive operator val- 
ued measure (POVM) A4 = {H x } having k detection oper- 
ators H x > which sums up to the identity J2 X H* = I- 
This measurement gives the correct answer x up to an error 
probability 

fe-i 
P = 1- X)p„ , R[IW* ) n)] := P[$\p(s,n),M] . (7) 

2 = 

Here we denote by P[$>\p(s, n),A4] the error probability in 
the readout of the cell $ given an input state p(s, n) and an 
output measurement M. . Now we are interested in minimizing 
this quantity over input and output. 

As a first step we fix the signal profile (s, n) and consider 
the minimization over input states and output measurements. 
This leads to the quantity 



P($|s,n) 



min P[$\p(s,n),M] 

p(s,n),M 



which is the minimum error probability achievable for a fixed 
signal profile (s,n). Note that there are some cases where 
the optimal output POVM is known. For instance if the out- 
put states p x (s, n) are pure and form a geometrically uniform 
, then the optimal detection is the square root mea- 




energy per cell, i.e., 

P($|n) = infP($|s,n) 



(9) 



Thus, given a memory with cell $, the determination of 
P($|n) provides the "optimal" readout of the cell at fixed en- 
ergy n. It is worth stressing that the minimization over the 
number of signals s is not trivial due to the constraint that 
we impose on the mean total energy (if instead of such re- 
striction one imposes a bound on the mean energy per signal, 
then the infimum is always achieved in the asymptotic limit 
of s — > oo). Also notice that we have put the word "opti- 
mal" between apostrophes, since the optimality of Eq. (0 is 
still partial, i.e., not including all the possible readout strate- 
gies. In fact, as we discuss in the following subsection, Bob 
can also consider the help of ancillary systems while keeping 
equal to n the mean total number of photons irradiated over 
the cell. 



B. Assisted readout of optical memories 

The optimality of Eq. (O is true only in the "unassisted 
case" where all the input modes are sent through the target 
cell. More generally, Bob can exploit an interferometric-like 
setup by introducing an ancillary "reference" system which 
bypasses the cell and assists the output measurement as de- 
picted in Fig. |2] In the "assisted case" we consider an input 
state p € T>(H® S ® / rl% r ) which describes s signal modes 
(Hilbert space T-if s ) plus a reference bosonic system with r 
modes (Hilbert space H% r ) 113011 . As before, the minimal en- 
ergy constraint corresponds to fixing the mean total number of 
photons irradiated over the target cell, i.e., n =Tr(phs) where 
hs is the total number operator acting over 'Hf ls Bill . We de- 
note by p = p(s,r,n) such a state, where we make explicit 
the number of signal modes s, the number of reference modes 
r, and the mean total number of photons n irradiated over the 
cell. Following the language of Ref. luOll . we also refer to 
p(s, r, n) as to a transmitter with s signals, r references, and 
signalling n photons B2I1 . 

Now, given a transmitter p(s,r,n) at the input of a target 
cell $ = {cj) x ,p x }, we have the output state 



Px(s,'< 



I® r )p(s,r,n) , 



(10) 



where the channel <fi x acts on each signal mode, while the 
identity / acts on each reference mode. This state is then 
measured by a POVM A4 = {H x } where n^ acts on the 
/g\ whole state space £>(%f s <g) H% r ). The error probability 
P[$>\p(s,r, n),A4] has the form of Eq. (O where now both 
state and measurement are dilated to the reference system. 
Thus, given a memory with cell $, the minimum error proba- 
bility at fixed signal energy n is given by 



As a final step, we keep the energy n fixed and we minimize 
over s, thus defining the minimum error probability at fixed 



P($|n)=inn min P[$\p(s,r,n),M] } , (11) 

s,r [p(s,r,n),M J 

where the minimization includes the reference system too. In 
general, we always consider the assisted scheme and the cor- 



responding error probability of Eq. ([TO . This clearly repre- 
sents a superior strategy for the possibility of using entan- 
glement between signal and reference systems. Clearly, the 
unassisted strategy is achieved back by setting r = and 

p{s,Q,n) = p(s,n). 
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FIG. 2: Assisted readout of an optical memory. Alice stores data 
in the cell by encoding a single-mode bosonic channel <f> x picked 
from the ensemble $. In general, Bob queries the cell by using a 
transmitter p(s, r, n) which describes s signal modes, irradiating n 
mean photons over the cell, plus r reference modes (bypassing the 
cell). The global output state p x (s, r, n) is detected by a quantum 
measurement M, which provides the correct answer x up to an error 
probability P[$|p(s, r, n), M]. 



III. THE SIMPLEST CASE: OPTICAL MEMORY WITH 
BINARY CELLS 

In general the solution of Eq. (TTTb is extremely difficult. 
In order to investigate the problem, the simplest possible sce- 
nario corresponds to an optical memory whose cell encodes 
two bosonic channels (binary cell) 12011 . The situation is par- 
ticularly advantageous when the channels are pure-loss and 
they are chosen with the same probability. This means to con- 
sider the binary channel ensemble 



$ = {4>u, Pu}u=0,l = {(t>0,P0,4>l-,Pl} I 



(12) 



where po — p\ = 1/2, and <f> u represents a pure-loss channel 
with transmission < k u < 1. In the Heisenberg picture, the 
action of <fi u on each signal mode is given by the map 



as 



u as - \/T" 



KuO.E , 



(13) 



where as is the annihilation operator of the signal mode and 
cie is the one of an environmental mode which is prepared 
in the vacuum state. For simplicity we can also denote the 
ensemble by using the short notation 



$ = {k ,Ki} 



(14) 



When the optical memory is read in reflection (which is usu- 
ally the case), then the two parameters kq and n\ represent 
the two possible reflectivities of the cell (so that unit reflec- 
tivity corresponds to perfect transmission of the signal from 
transmitter to receiver). 



Given a transmitter p(s,r,n) at the input of the binary 
cell $, we have two equiprobable outputs, p Q (s,r,n) and 
pi (s, r, n). In this case the optimal measurement corresponds 
to the projection onto the positive part of the Helstrom matrix 
Po{s, r, n) — pi(s, r, n) fljj] . As a result, the error probability 
for reading the binary cell $ using the transmitter p(s, r, n) is 
given by 

P[$\p(s,r,n)] = - jl- -D[p (s,r,ri), pi(s,r,n)] 

(15) 
where D is the trace distance [1]. This expression has to be 
optimized on the input only, so that we can write 



P(®\n)=m£i min P[$|p(s, r, n)] 

s,r I p(s,r,n) 



(16) 



which is the minimum error probability at fixed signal energy. 
This quantity clearly provides the maximum information per 
cell at fixed signal energy, which is given by 



where 



I($\n) = l-H[P($\n)], (17) 



H(x) = -xlog 2 x-(l-x)log 2 (l-x) (18) 



is the binary formula of the Shannon entropy. 

Even in this simple binary case the solution of Eq. ( fT6l ) 
is very difficult. However we can provide remarkable lower 
bounds if we restrict the minimization to some suitable class 
of transmitters. An important class is the one of the classical 
transmitters, since they encompass all the optical resources 
used for the readout of optical memories in today's storage 
technology. Furthermore, this class can be easily character- 
ized. Given a transmitter p(s, r, n), we can write its Glauber- 
Sudarshan representation 13311 



p(s, r,n)= / d 2s a d 2r (3 P(a, /3) a(a) g> >y(J3) , (19) 

where a = (ai, • • • , a s ) T and /3 — (/3i, • • • , f3 r ) T are vec- 
tors of complex amplitudes, 

S T 

a(a) = (g) \oi)s{ai\ , l(fi) = (g) |A)r(A| , (20) 
i=l i=l 

are multimode coherent states, and the P-function P(a, /?) is 
a quasi-distribution, i.e., normalized to one but generally non 
positive 113311 . In terms of the P-function, the signal energy 
constraint reads 



Jd 2s ad 2r pP(a,/3)J2 



cti\ = n 



(21) 



Now we say that p(s, r, n) is classical (nonclassical) if the P- 
function is positive (non positive). Thus if the transmitter is 
classical, denoted by p c (s, r, n), then it can be represented as 



a probabilistic mixture of coherent states. The simplest ex- 
amples of classical transmitters are the coherent state trans- 
mitters, that we denote by p co h(s, r, n). These are defined by 
singular P-functions 

P{a,(3) = 6 2s {a-a)6 2r {/3-P), (22) 

so that they have the simple form 

pcoh{s, r, n) = a(a) <g> j(J3) . 

Examples of nonclassical transmitters are constructed using 
squeezed states, entangled states and number states [24, 25]. 

As shown in Ref. 12011 . by restricting the optimization to 
classical transmitters, we can compute the upper bound 



P(*|n) < P c (*|n) 



inf < min P[<&\p c (s,r,n)] 

8,r I p c (s,r,n) 



which is given by 

Pc(®\n) = - 



(23) 



(24) 



This bound can be reached by a coherent state transmitter 
Pcoh(^,0,n) = | y/n)s(y/n\, i.e., a single-mode coherent state 
with mean number of photons equal to n. The error probabil- 
ity P c (<t\n) of Eq. d24i i or, equivalently, the mutual informa- 
tion 



I c ($|n) = 1 - H[P c (?\n)] 



(25) 



is known as "classical discrimination bound". 

Alternative (and better) bounds can be derived by resorting 
to nonclassical transmitters 12011 . As a prototype of nonclassi- 
cal transmitter we consider the EPR transmitter IU8I1 . which is 
composed by s pairs of signals and references, entangled via 
two-mode squeezing. This transmitter has the form 



p epr {s,s,n) = |£) (£;f 



(26) 



where |£) (£| is a two mode squeezed vacuum (TMSV) state, 
entangling one signal mode S with one reference mode R. In 
the number-ket representation, we have 12411 



|0 = (coshO" 1 X>anliO f 



IR ' 



(27) 



where the squeezing parameter £ quantifies the signal- 
reference entanglement and gives the energy of the signal by 
sinh 2 £, Since this transmitter involves s copies of this state, 
we have to impose 



£ = arc sinh \j — , 



(28) 



in order to have an average of n total photons irradiated 
over the cell. Given an EPR transmitter p epr (s, s, n) at the 
input of the binary cell $, we have an error probability 



P[$\p epr (s, s, n)]. By optimizing over the number of copies 
s, we define the upper bound 

P(*|n) < Pe P r($\n) :=MP[$>\p epr (s,s,n)} . (29) 

This bound represents the maximum information which can 
be read from the binary cell $ by using an EPR transmitter 
which signal n mean photons. This quantity can be estimated 
using the quantum Battacharyya bound and its Gaussian for- 
mula |4J] . After some algebra, we get J20l 12611 



P epr ($|n) < B 



exp(— um) 



(30) 



where 



u := 2^/kqKi - y/(l - «o)(l - «i) • (31) 



A. Quantum versus classical reading 



Because of the potential implications in information tech- 
nology, it is important to compare the performances of clas- 
sical and nonclassical transmitters. The basic question to ask 
is the following [|20ll : for fixed signal energy n irradiated over 
a binary cell $, can we find some EPR transmitter able to 
outperform any classical transmitter? In other words, this is 
equivalent to show that P epr ($|n) < P c ($|n) and a sufficient 
condition corresponds to prove that B < P c ($|n). Thus, by 
using Eqs. (T24l and (l30l , we find that for signal energies 



n > n t h 



2 In 2 



2 - K0 — Kl - 2y/(l - K )(l - «l) 



(32) 



it is always possible to beat classical transmitters by using an 
EPR transmitter 112011. For high reflectivity k± ~ 1 and kq < 
K\ the threshold energy nth can be very low. In the case of 
"ideal memories", defined by kq < Ki — 1, the bound of 
Eq. ( f30b can be improved. In fact, we can write 



P epr (*|n) < 



exp[-2n(l - y/iio)] 



(33) 



and the threshold energy becomes nth = 1/2 02011 . Thus, 
for optical memories with high reflectivities and signal ener- 
gies n > 1/2, there always exists a nonclassical transmit- 
ter able to beat any classical transmitter. In the few-photon 
regime, roughly given by 1/2 < n < 10 2 , the advantages 
of quantum reading can be numerically remarkable, up to one 
bit per cell. The implications have been thoroughly discussed 
in Refs. J20l 12611 . It is important to say that these advan- 
tages are also preserved if thermal noise is added to the ba- 
sic model. This noise can describe the effect of stray photons 
hitting the memory from the background and other decoher- 
ence processes occurring in the reading device. Formally, this 
means to extend the problem from the discrimination of pure- 
loss channels to the discrimination of more general Gaussian 
channels 12611. 



A supplementary analysis of quantum reading has also 
shown that its advantages persist if we consider more ad- 
vanced designs of memories where information is written 
on and read from block of cells (multi-cell/block encoding). 
Block encoding allows Alice to introduce error correcting 
codes which make Bob's readout flawless up to some meta- 
data overhead. By resorting to the Hamming bound and the 
Gilbert- Varshamov bound, Ref. 12611 showed that EPR trans- 
mitters enable the low-energy flawless readout of classical 
memories up to a negligible error correction overhead, con- 
trarily to what happens by employing classical transmitters. 
In the following section, we develop the idea of block encod- 
ing in the most general scenario, i.e., for arbitrary classical 
memories. Then, by sending the size of the block to infinite 
(Sec. |V), we will be able to introduce the notion of quantum 
reading capacity of a classical memory. 



IV. GENERAL MODEL OF MEMORY: MULTI-CELL 
ENCODING 

The writing of a memory is based on channel encoding 
which generally may involve a block of m cells. A first triv- 
ial kind of block encoding is just based on independent and 
identical extractions. As usual, Alice encodes a fc-ary vari- 
able X = {x,p x } into an ensemble of quantum channels 
<& = {4> x ,p x }. Then, she performs m independent extrac- 
tions from X, generating an m-letter sequence 



(xi 



.) 



(34) 



with probability p x = p Xl ■ ■ -p Xm . This classical sequence 
identifies a corresponding "channel-sequence" 

(f>x ■■= 0xi <8> ■ ■ • <8> <t>x m , (35) 

which is stored in the block of m cells. 

In a more general approach, Alice adopts a classical code. 
This means that Alice disposes a set of ?7i-letter codewords 
{x°, • • • , x 4 , • • • , x' -1 } with I < k m . A given codeword 

K i = (x\,--- ,o4), (36) 



is chosen with some probability p x 
sponding "channel-codeword" 



and identifies a corre- 



(37) 



Thus, in general, Alice encodes information in a block of m 
cells by storing a channel-codeword, which is randomly cho- 
sen from the ensemble {</> x ; , p x i } where % = 0, • • • ,1—1. 

The most general strategy of readout can be described as 
a problem of "parallel discrimination of quantum channels", 
where Bob probes the entire block in a parallel fashion and 
detects the output via a collective quantum measurement. In 
order to query the block, Bob uses s signal systems per cell 
besides other supplemental r reference systems for the benefit 
of the output measurement. The whole set of ms + r systems 
is described by an arbitrary multipartite state p (see Fig. [3]). 
At the output of the block, Bob has 
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FIG. 3: Memory model with block encoding. In order to write 
data, Alice encodes a channel-codeword cj) x i in a block of m cells. 
To read the data, Bob uses a suitable transmitter and receiver solv- 
ing a problem of parallel channel discrimination. The transmitter is 
an arbitrary multipartite state p which probes the entire block by in- 
putting $ systems per cell plus sending additional r systems directly 
to the receiver. The output state p^t is detected by an optimal collec- 
tive measurement which provides the correct answer x ! up to some 
error probability P err . In the uncostrained readout, P err goes to 
zero and Bob retrieves all the information l/ mal from the block. If 
the readout is constrained, as in the case of optical memories at fixed 
signal energy n, then P err is nonzero. In this case, Bob retrieves a 
fraction of the information I < J/ m ax or, equivalently, he retrieves 
all the information if Alice suitably increases the size of the block 
while keeping -ff max as constant. 



where the identity acts on the reference systems, while 



(39) 



acts on the signal systems. This state is detected by a collec- 
tive quantum measurement, i.e., a general POVM with I de- 
tection operators with outcome i corresponding to codeword 
x\ Clearly, the main goal for Bob is to optimize both input 
state and output measurement in order to retrieve the maximal 
information from the block. 

It is intuitive to understand that, without constraints, Bob is 
always able to retrieve all the information from the block, i.e., 



/-i 



H t] 



= - > ,Px>l0gP x 



=0 



(40) 



•I®1(P), 



(38) 



bits of information. However this is no longer the case if we 
impose constraints on Bob's physical resources. As we know, 
when we consider optical memories (bosonic setting), the op- 
timization must be constrained in the input energy, in partic- 
ular, by fixing the mean total number of photons n irradiated 
over each cell. In this case, if we consider low values of n, the 
measurement will be affected by non-negligible error proba- 
bility P err and the information retrieved will be some value 
I = I(n) between and H max . In other words, data will be 
read with an average rate of R(n) = ?ri~ 1 /(n) bits per cell. 

It is important to note that, in the block-encoding model, an 
equivalent approach consists of making the readout flawless 



by increasing the error correction overhead in the block. In 
other words, for a given signal energy n, we can determine the 
minimal size m = m(n) of the block (and the corresponding 
optimal classical code) which makes the error probability P err 
negligible (i.e., reasonably close to zero 112610 . In this case, 
the readout is flawless and the block provides all the H msx 
bits of information. As a result, the rate now takes the form 
R(n) = Mn)]- 1 ^^. 

It is clear that, given an arbitrary block encoding {cf> x i , p x i } 
and an arbitrary multipartite transmitter which irradiates n 
mean photon per cell (see Fig. |3), the computation of the rate 
R(n) is extremely difficult. However we can face the problem 
if we consider transmitters which are separable with respect to 
the different cells (more exactly, in tensor product form) and 
taking the limit for infinite block (m — > oo). This allows 
us to introduce a simple description of the memory (similar 
to the single-cell scenario) and, most importantly, to use the 
Holevo bound as quantifier for the readable information, i.e., 
as asymptotic rate R(n). Then, the optimization of this rate 
over the transmitters enables us to define the quantum reading 
capacity of the memory. 
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FIG. 4: Limit for infinite block. A memory can be described as a 
large (approximately infinite) encoding block, where each cell en- 
codes a marginal ensemble $ = {<f> x ,p x }. In order to read the 
memory, Bob uses an multi-copy transmitter p(s, r)®°° = p(s, r) ® 
p(s, r) ® ■ ■ • , where each copy p(s, r) probes a different cell using 
s signals generally coupled with other r references. All the outputs 
from the block are collectively detected by an optimal quantum mea- 
surement which reconstructs the asymptotic channel codeword. 



V. LIMIT FOR INFINITE BLOCK: QUANTUM READING 
CAPACITY 

Digital memories typically store a great amount of data. 
This means that an average memory is composed of a large 
number of encoding blocks. In principle, we can also describe 
the memory as a single large block of cells where Alice stores 
data by encoding a very long channel-codeword </> XI chosen 
with some probability p x i . Considering the whole memory 
as a large encoding block allows us to re-introduce a single- 
cell description. In fact, in the limit for m — > oo, each cell 
can be described (on overage) by a marginal ensemble of 
quantum channels $ = {cj) x ,p x } encoding a corresponding 
marginal variable X = {x,p x }. Thus, independently from 
the actual classical code used to store information, the de- 
scription of a large classical memory can always be reduced 
to its marginal cell, corresponding to a marginal ensemble of 
channels $ = {4> x ,p x }. 

Despite this asymptotic simplification, the readout process 
is still too difficult to be treated if we consider arbitrary mul- 
tipartite states, i.e., generally entangled among different cells. 
Thus we restrict the readout to input states which are tensor 
products. This means that Bob inputs an oo-copy state 



where 



p(s,r) 



p(s,r)®p(s,r) 



(41) 



where the single-copy p(s,r) £ X>(/Hf s <g> U% r ) describes 
s signal systems sent through a target cell plus additional r 
reference systems (see Fig. |4j. Given the oo-copy transmit- 
ter p(s,r)^' co at the input of a memory with marginal cell 
"5 = {4>x,Px}, the output is still in a tensor product form. 
The average output of each cell is described by a marginal 
ensemble of states 



{Px(s,r),p x } , 



p x (s,r) = («^ s <g> I® r )[p(s,r)} 



(43) 



By applying an optimal collective measurement on all the out- 
puts, the maximum information per cell that can be retrieved 
is given by the Holevo bound 



x{£) = s 



^2p x p x {s,r) 



-J2PxS[ Px (s,r)} (44) 



where S denotes the von Neumann entropy. The achievability 
of x(£) is assured by the Holevo-Schumacher- Westmoreland 
theorem Jg]. Here it is important to note that the asymp- 
totic readout can be flawless. In other words, for a given 
marginal $ = {(j> x ,Px} there is always an asymptotic block 
code {4> x i,p x i}, giving that marginal, which allows the re- 
ceiver to give the correct answer x' with asymptotically zero 
error. 

Given a memory with marginal cell $, the Holevo infor- 
mation of Eq. ( f44l > depends on the input state p(s, r) only. 
This means that it can be represented as a conditional Holevo 
information, that we denote by x[3>|/?(s, r)]. In other words 
X[$|p(s,r)] represents the maximum information per cell 
which can be read from a memory with marginal cell $ if 
we use the transmitter p(s, r). The crucial task here is the op- 
timization of x[$|/?(s, r)] over the transmitter. As a first step, 
we can consider the readout capacity for fixed number of input 
systems s and r, i.e., 



C($\s,r) 



rnaxx[$|p(s,r)] , 



(45) 



Now, by optimizing Eq. d45l ) over the number of input sys- 
tems, we can define the unconstrained quantum reading ca- 
pacity of the memory [350 



(42) 



C(<&) = sup max x[<I>|p(s, r)] 

s.r p(s,r) 



(46) 



Since it is unconstrained, this capacity can be greatly sim- 
plified and trivially computed. First of all the maximization 
can be reduced to pure transmitters ijj(s, r) as a simple conse- 
quence of the convexity of the Holevo information B34fl (see 
Appendix I A II for more details). Then, also the use of the ref- 
erence systems can be avoided. In other words, it is sufficient 
to consider the unassisted capacity where we maximize over 
i/j(s, 0) = ip(s). Furthermore the supremum is achieved in the 
limit for s — ► +00, i.e., we can write 

C($) = lim maxx[$|V'0)] • ( 47 ) 

This quantity is the maximal possible since it equals the 
amount of information stored in the marginal cell of the mem- 
ory. This is given by the Shannon entropy of the marginal 
variable X = {x,p x } encoded by the marginal ensemble 
*£ = {4>x,Px}- In other words, we have 



C($) = H(X) = -J2p^°SP* 



(48) 



The proof is trivial (see Appendix IA 2l for details). 

The notion of quantum reading capacity is non-trivial only 
in the presence of physical constraints. This is what happens 
in the bosonic setting, where optical memories are read by 
fixing the input signal energy. Thus, let us consider an optical 
memory with marginal cell $ = {4> x ,Px} where 4> x repre- 
sents a single-mode bosonic channel. As transmitter, now we 
consider an oo-copy state 



p(s, r, n) 



p(a,r,n) <g> p(s,r,n) 



(49) 



where p(s, r, n) £ T>(l-L® s ® T~L% r ) describes s signal modes, 
irradiating n mean photons on a target cell, plus additional r 
reference modes bypassing the cell. At the output we have an 
infinite tensor product of states of the form 



Px(s,r,n) 



.I® r )Ws,r,n)}, 



(50) 



which are detected by an optimal collective measurement. In 
this way, Bob is able to retrieve an average of %[$|p(s, r, n)] 
bits per cell. Now, we must optimize this quantity over the 
input transmitters by taking the signal energy n fixed. This 
constrained optimization leads to the definition of the quan- 
tum reading capacity of the optical memory 

C($|n) = sup max x[®\p( s , r , n )] ■ (51) 

s.r P{s,r,n) 

This capacity represents the maximum information per cell 
which is readable from an optical memory $ by irradiating 
n mean photons per cell. The computation of Eq. ( BTT i is not 
easy at all. As a matter of fact we are only able to provide 
lower bounds by restricting the class of transmitters involved 
in the maximization. We do not even know if the optimal 
transmitters are pure or mixed. 

Let us consider a set (or "class") V of pure transmitters 
ip(s,r,n) which are characterized by some general property 
which does not depend on s, r and n (for instance, they could 
be constructed using states of particular kind, such as coherent 



states). Then we can always construct the mixed-state trans- 
mitter 



(52) 



p(s, r,n)= / dy p y i/) y (s, r, n) , 



where p y ^ 0, J dy p y = 1, and ^p y (s,r,n) e V. Clearly 
the set of mixed-state transmitters identifies a larger class A 
which includes V. Now, we can define a lower-bound to 
C($|n) by optimizing over the class A, i.e., 

C($|n) > CU($|n) = sup max x[^\p(s,r,n)} . (53) 

s.r p(s,r,n)t=A 

Similarly we can consider the further lower-bound 

CU($|n) > C-p($|n) = sup max x[$|V->(s, r, n)] . 

S] r 4>(s,r,n)eV 

(54) 
Here we first ask: is there some class V that allows to put an 
equality in Eq. (O, i- e -, Op($|n) = CU($|n)? Then, is it 
possible to extend this class to all the pure transmitters, so that 
CV>($H=C($|n)? 

Unfortunately we are not able to answer the second ques- 
tion, so that the issue of the purity of the optimal transmitters 
remains unsolved. However, we are able to find classes for 
which Cp($|n) = CU($|n). For this sake, a sufficient crite- 
rion is the concavity of Gp($|n). 

Lemma 1 IfC-p(<&\n) is concave in n, then we have 

Cp(*|n) = CU(*|n) . (55) 

Proof. Let us consider the transmitter of Eq. d52l , whose sig- 
nal energy (mean number of photons) can be written as 



(56) 



= I dy p y n y , n y = (ip y \h\ip y } 



Given this transmitter at the input of a marginal cell $, we can 
bound the conditional Holevo information 



X[®\p(s,r,n)} ^ / dypyx[m> v ] ( 57 ) 

< fdy Py C v (^\n y ) (58) 



< C v U\ J dy p y ny) (59) 

= C v ($|n) , 



(60) 



where we have used the convexity of % in the first inequal- 
ity d57l ), the definition of C-p($|n) in the second inequal- 
ity J58l and its concavity in the last inequality J59l ). It is clear 
that Eqs. d57ll-(l60ll hold for every p(s, r,n) e A and every s 
and r. As a result, we can write 

sup max x[®\p( s , r , n )} = CU($W < Cv ($|n) , 

s.r p(s,r.n)eA 

(61) 

which, combined with Eq. d54l ), gives the result of Eq. ( 1331 ). ■ 

In the following section we show that an important class V 

for which C-p ($|n) is concave is the one of the coherent-state 
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transmitters. This means that Cp($|n) = CU($|n), where A 
is the class of the classical transmitters (constructed by con- 
vex combination via the P-function). Thanks to this result we 
can compute an analytical bound for the readout performance 
of all the classical transmitters, that we call "classical read- 
ing capacity". This capacity represents the multi-cell gener- 
alization of the classical discrimination bound of Sec. Hill and 
provides a simple lower bound to the quantum reading capac- 
ity. In Sec.|Vl]we compute its analytical formula for the most 
basic optical memories. Then, as we will show in Sec. |VH| 
this classical bound can be easily outperformed by nonclassi- 
cal transmitters, thus proving its separation from the quantum 
reading capacity. 



VI. CLASSICAL READING CAPACITY 

Let us consider an optical memory which is the multi-cell 
generalization of the binary model described in Sec. [HI] In 
the single-cell model of Sec. [TTTJ information was written in 
each cell in an independent fashion, by encoding one of two 
possible pure-loss channels, <fro and (\>\ (binary cell). Here we 
consider the multi-cell version, where Alice stores a channel 
codeword in the whole optical memory regarded as an infi- 
nite block. In particular, the block encoding is such that the 
marginal cell is described by a binary ensemble 



where H is the binary Shannon entropy and 



* = {<t>o,p,<t>\A -p} 



(62) 



where < p < 1 and <f> u is a pure-loss channel with transmis- 
sion k u . Alternatively, we can use the notation 

$ = {no,p,Ki,l-p} . (63) 

Given this kind of memory, we consider the input 

p c (s, r, n)®°° = p c (s, r, n) ® p c (s, r,n)® ■■ ■ , (64) 

where p c (s, r, n) is an arbitrary classical transmitter with s 
signals, r references and n mean photons. The average in- 
formation which can be read from each cell is provided by the 
Holevo quantity x[$\Pc(s, r, n)]. Now, by optimizing over the 
classical transmitters we can define the lower-bound 

C($|n) > C c ($|n) = sup max x[$\pc{s,r,ri)] , (65) 

s.r p c (s,r,n) 

which represents the classical reading capacity of the optical 
memory $. This capacity represents the multi-cell version of 
the classical discrimination bound of Sec. [HI] As before, we 
can provide a simple analytical result. 

Theorem 1 Let us consider an optical memory with binary 
marginal cell $ = {no,p, K\, 1 — p} which is read by a clas- 
sical transmitter signalling n mean photons. Then, the maxi- 
mum information per cell which can be read is asymptotically 
equal to 



1 1 



£ = 2 + 2V 1 " 4p(1 ~ p) t 1 " e ~ n{VTri ~ 



(67) 



In particular, the bound C c ($\n) can be reached by using a 
coherent-state transmitter p co ^(l,0,rt) = \^/n)s(\/n\, i.e., a 
single-mode coherent state with n mean photons. 

Proof. Let us consider the class V = coh of coherent-state 
transmitters p co h(s,r,n). By convex combination we con- 
struct the class A = c of the classical transmitters p c (s, r, n). 
The first step of the proof is the computation of C co h(^\n), 
i.e., the readout capacity restricted to coherent state transmit- 
ters. We first prove that C coh {§ \ n) = x[$\Pcoh(l,0,n)], 
i.e., the optimal coherent-state transmitter is the single-mode 
coherent state \y/n)s{\fn\. Then, we analytically compute 
xi^lPcohO-, 0, n)]. Since this quantity turns out to be concave 
in n, we can use LemmaQ]and state C co h(&\n) = C c (Q\n), 
thus achieving the result of the theorem. 
Given a coherent-state transmitter 

p coh (s,r,n) = a{a) <8>7(/3) 

s r 

= ®|ai}s<Oi|®(g)|/8i)H<AI (68) 

i=l i=\ 

at the input of the cell $, we have the output 

Pu = <ff[o{a)\® 1 {P) 

s r 

= ®M|oi>s(ail)®®IAMAI 

i=l i=l 

s r 

= (g) \V^^)s(V^^\ ® lAMftl , (69) 

i=\ t=l 

which is still a multimode coherent state. This is a simple 
consequence of the fact that <fio and <f>\ are pure-loss channels. 
Since we are computing the Holevo information on the output 
ensemble, we have the freedom to apply a unitary transforma- 
tion over p u . By using a suitable sequence of beam splitters 
and phase-shifters we can always transform p u into the state 



i>r+s— 1 



(70) 



Then, since the Holevo information does not change under the 
adding of systems, we can trace the r + s — 1 vacua and just 
consider the single-mode output state 



Pu = WK u n)s(y/K u n\ 



(71) 



This can be achieved by considering a single-mode coherent 
state transmitter 



p coh (l,0,n) = \Vn)s(Vn\ 



(72) 



C c (*|n) = .fftf) , 



(66) 



at the input of the pure-loss channel <p u . For fixed marginal 
cell $ and fixed input energy n, the reduction from the multi- 
mode input of Eq. d68l l to the single-mode output of Eq. d72l ) 
is always possible, independently from the actual number of 
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systems, s and r, and the specific form of the transmitter 
Pcoh{ s , r i n ). Thus, we can write 



C coh {®\n) = sup max x[$\pcoh(s,r,ri)} 

s,r Pcoh{s,r,n) 

= x[*|peoh(l,0,n)]. 



(73) 



In other words, the optimal coherent-state transmitter is the 
single-mode coherent state \y/n)s(\/n\. The next step is the 
analytical computation of x[&\Pcoh(l, 0, ri)]. After some Al- 
gebra we get 

xl$\pcoh(l,0,n)} = H(0, 

where H is the binary formula of the Shannon entropy and 
£ = £,( K o, K i,Pi ri) is given in Eq. ( f67b . One can easily check 
that H(£) is a concave function of n, for any kq, K\, and p. 
Since C co ^($|n) = i?(£) is concave in the energy n, we can 
apply Lemma [TJ by setting V = coh and .A = c. Thus we 
get C c (l>|n) = C' co / 1 (l'|?i) = i?(£) which is the result of 
Eq. d66i l. It is clear that the optimal classical transmitter co- 
incides with the optimal coherent-state transmitter which is 
given by p coh {l,0,ri). ■ 

It is interesting to compare the single-cell and multi-cell 
classical discrimination bounds, in order to estimate the gain 
which is provided by the parallel readout of the cells. For 
a direct comparison, let us set p = 1/2, so that the binary 
cell $ is described by $ = {«o, 1/2, «i, 1/2} = {k ,ki}. 
Then, we compare the maximum information which achiev- 
able by using classical transmitters in the multi-cell readout, 
i.e., the classical reading capacity C c ($|n), with the maxi- 
mum information which is achievable by classical transmit- 
ters in the single-cell readout, i.e., the classical discrimination 
bound I c (<&\n) given in Eq. d25l l. As shown in Fig.|5]the ad- 
vantage is quite evident. 



Q. 




FIG. 5: Maximum number of bits per cell read by classical trans- 
mitters as a function of the signal energy n (mean number of 
photons). We compare the two classical discrimination bounds: 
C c (&\n) (multi-cell readout, solid line) and 7 c ($|n) (single-cell 
readout, dashed line). We consider a memory with binary marginal 
cell <E> = {acq, «i} where kq = 0.5 and «i = 0.9. 



In the following Sec. IVIII we will construct examples of 
nonclassical transmitters which are able to outperform the 
classical reading capacity. This will prove the separation be- 
tween the quantum reading and the classical reading capaci- 
ties, thus showing the advantages of quantum reading in the 
multi-cell scenario. 



VII. NONCLASSICAL TRANSMITTERS 

As before, let us consider an optical memory with binary 
marginal cell $ = {k , K\\. This time we assume that it is 
read by using a nonclassical transmitter p nc (s, r, n). Since we 
are in the asymptotic multi-cell scenario, we clearly assume an 
oo-copy input p nc (s, r, n)<Eip nc (s, r,n)<E>- ■ ■ together with an 
optimal collective measurement of the output. The maximum 
number of readable bits per cell is given by the conditional 
Holevo information x[$|p nc (s, r, ri)]. Now we ask: is this 
quantity bigger than the classical reading capacity C c ($|n)? 

The first design of nonclassical transmitter is the EPR trans- 
mitter p epr (s, s, ri) = |£) (£| s which has been first discussed 
in Sec. [TTTJ In order to beat classical transmitters, it is suffi- 
cient to consider p epr (l, l,ri) — |£) (£|, i.e., a single TMSV 
state per cell. This means that we have one signal mode S, 
irradiating n mean photons over a target cell, which is entan- 
gled with one reference mode R. To quantify the advantage 
we consider the information gain 

G = X [*|/»epr(l ) l,n)]-C c ($|n) 

and check its positivity. If G > then the EPR transmit- 
ter p epr (l, 1, n) beats all the classical transmitters, retrieving 
G bits per cell more than any classical strategy. As shown 
in Fig. [6] we have G > in the regime of low photons and 
high reflectivities (i.e., kq or K\ close to 1). This is the typ- 
ical regime where the quantum reading of optical memories 
is advantageous, as also investigated in the single-cell sce- 
nario 12011. 
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FIG. 6: (Color online) Information gain G versus reflectivities, ko 
and Ki, for n = 5 (left panel) and n — 1 (right panel). Here G pro- 
vides the number of bits per cell which are gained by the single-copy 
EPR transmitter |£) (£| over all the classical transmitters in the read- 
out of an optical memory with marginal cell $ = {^o, fti}. Note 
that the highest values of G occur for Ko or ki close to 1 (high re- 
flectivities). 

As evident from Fig. [6] the best situation corresponds to 
having one of the two reflectivities equal to 1, i.e., for an "ideal 
memory" $ = {kq < K\, «i = 1}. Given such a memory, we 
explicitly compare the information read by an EPR transmitter 
Xepr — x[^|Pepr(l, 1, Ti)] with the classical reading capacity 
C c ($|n) at low signal energy. As shown in Fig.|7]for n = 1, 
the EPR transmitter is always able to beat the classical bound. 

It is important to note that we can construct other simple ex- 
amples of nonclassical transmitters which can outperform the 
classical reading capacity. An alternative example of nonclas- 
sical transmitter can be taken again of the form p nc (l, 1, ri) 
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FIG. 7: Number of bits per cell as a function of Ko, for «i = 1 
(ideal memories) and n — 1 mean photons per cell. We compare 
the classical reading capacity (dotted line) with the Holevo informa- 
tion retrieved by various nonclassical transmitters: EPR transmitter 
(dash-dotted line), NOON state transmitter (dashed line) and Fock 
state transmitter (solid line). 



and corresponds to the NOON state H361I37H 

\NOON) = 2- 1 / 2 (|2n) s |0) i? + \0)s\2n) R ) , (74) 

where signal and reference are again entangled. A further ex- 
ample of nonclassical transmitter is of the form p nc (l, 0, n), 
i.e., not involving the reference mode. This is the Fock state 



\n) s = (n\)-^(alr\0) i 



(75) 



As shown in Fig. [7] these transmitters can beat not only the 
classical reading capacity but also the EPR transmitter |£) (£| 
for low values of kq. Recently, these kinds of transmitters 
have been also studied by Ref. 112111 in the basic context of 
quantum reading with single-cell readout. 

It is interesting to compare the performances of all these 
transmitters in the low-energy readout of optical memories 
with very close reflectivities. This is shown in Fig. [8] for 
^1 — Ko = 0.01 and n — 1. The EPR transmitter |£) (£| is 
optimal almost everywhere, while the classical bound beats 
the other nonclassical transmitters for low values of «i (and 
Ko). This is also compatible with the result of optimality of the 
TMSV state for the problem of estimating the unknown loss 
parameter of a bosonic channel 0811 . As evident from Fig. [8] 
the bigger separation from the classical bound occurs for high 
reflectivities, i.e., K\ close to 1. 

It is also interesting to see what happens in the regime of 
low reflectivity by considering a binary marginal cell $ = 
{ko, «i} with «o = 0. For this comparison we introduce an- 
other nonclassical transmitter of the form p nc (l, 0, n). This 
is the squeezed coherent state \a,£) = D(a)S(£)\0), where 
D(a) is the displacement operator and S(£) the squeezing op- 
erator J24I1 . The squeezed coherent state is chosen with the 
squeezing orthogonal to the displacement direction. Then we 
choose two real parameters, a and £, which are optimized un- 
der the condition a 2 + sinh £ = n, imposed by the mean 
photon-number constraint. As shown in Fig. |9j the presence 
of squeezing is sufficient to outperform the classical reading 
capacity in the regime of low reflectivity. However, better per- 
formances can be achieved by the Fock state considering high 
values of «i. 
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FIG. 8: Number of bits per cell as a function of ki, for Ki — kq = 
0.01 (close reflectivities) and n — 1. We compare the classical read- 
ing capacity (dotted line) with different nonclassical transmitters: 
EPR transmitter (dash-dotted line), NOON state transmitter (dashed 
line) and Fock state transmitter (solid line). 
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FIG. 9: Number of bits per cell as a function of ki, for no = and 
n — 1. We compare the classical reading capacity (dotted line) with 
the squeezed coherent state transmitter (dashed line) and the Fock 
state transmitter (solid line). 



From the previous analysis it is evident that, in the regime 
of low photon number (down to one photon per cell), we can 
easily find nonclassical transmitters able to beat any classical 
transmitters, i.e., the classical reading capacity. This is par- 
ticularly evident for high reflectivities (ki or ko close to 1). 
Thus, for the most basic optical memories, classical and quan- 
tum reading capacities are separated at low energies. In other 
words, the advantages of quantum reading are fully extended 
from the single- to the optimal multi-cell scenario. 

At this point a series of important considerations are in or- 
der. First of all, note that we have only considered nonclas- 
sical transmitters irradiating one signal per mode (entangled 
or not with a single reference mode), i.e., transmitters of the 
kind p nc (l, 0, n) or p nc (l, 1, ri). The reason is because these 
transmitters are sufficient to beat the classical bound. How- 
ever, better performances can be reached by optimizing over 
the number signals and references. In the case of the EPR 
transmitters, we expect that p epr (2, 2, n), which is composed 
of two TMSV states signalling n/2 mean photons each, is 
able to outperform p epr (l, 1, n), i.e., a single TMSV state sig- 
nalling ?i mean photons. This is shown in Fig.[l0]for the case 
of an ideal memory and n = 1 mean photons. This advan- 
tage could further improve for EPR transmitters p epr (s, s, n) 
with higher values of s. For this reason, in order to reach the 
quantum reading capacity, it is necessary to optimize over an 
arbitrary number of signal and reference modes, as foreseen 
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by the general definition of Eq. ( f5Tb . 




FIG. 10: Number of bits per cell as a function of Ko, for «i = 1 
(ideal memory) and n — 1. We compare the classical reading capac- 
ity (dotted line) with two different EPR transmitters: p epr (l, l,n) 
(lower solid line) and p ep r-(2, 2, n) (upper solid line). 

Another important consideration is related to the practical 
realization of quantum reading. In order to be experimentally 
feasible, the detection scheme should be as simple as possi- 
ble. For this reason, it is interesting to compare the classical 
reading capacity (which refers to the general multi-cell read- 
out) with the performances of EPR transmitters in the single- 
cell scenario, where each cell is detected independently from 
the others. Thus, we consider an ideal memory $ = {kq < 
Ki,K\ = 1} which is irradiated by a few mean photons per 
cell (in particular, we can consider n = 5). Given this mem- 
ory, we compare the optimal performance C c ($|n) of classi- 
cal transmitters assuming the multi-cell readout (asymptotic 
collective measurement) with the performance of EPR trans- 
mitters p ep r(s, s, n) assuming the single-cell readout (individ- 
ual cell-by-cell measurements). The latter quantity is given by 
the mutual information 



l epr ($\s,n) = l-H{P[^\p epr {s,s,n)]}, 



(76) 



where H is the binary formula for the Shannon entropy and 
P[<b\p epr (s, s,n)} is the error probability of the single-cell 
readout. One can compute the upper bound 



P[<S>\ Pepr (s lSl n)}<e:=- 



1+- 1 

s 



, (77) 
which provides a lower bound for the mutual information 



J r epr (*|*,n)>Q(fc|a,n) := 1 - #(0) 



(78) 



Thus, Q(<&\s, n) provides the minimum number ofbits per cell 
which are read by an EPR transmitter p epr {s, s, n). For fixed 
signal energy n, it is trivial to check that this quantity is in- 
creasing in s. This means that for any integer s we have 



Q($|l,n) < Q($\s,n) < Q($|oo,n) 



(79) 



where Q($|l,?i) corresponds to a single energetic TMSV 
state p epr (l,l,n) while Q($|oo,n) corresponds to 
p epr (oo, oo, n), i.e., infinite copies of TMSV states with 
vanishing energy. The quantity Q($|oo,n) is computed by 
taking the limit for s — > oo in Eq. ( ITTl i. In this limit, we have 



— >• 0, with 6 given in Eq. d33l . In the upper panel of Fig.fTTI 
we explicitly compare the two extremal values Q($| 1, n) and 
<2($|oo, n) with the classical reading capacity C c ($|n). As 
we can see, the single-cell quantum reading is able to beat the 
asymptotic multi-cell classical reading. 
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FIG. 11: Upper panel. Number of bits per cell as a function of reo, 
for ki = 1 (ideal memory) and n — 5. We compare the classical 
reading capacity C c ($|n) (multi-cell readout, dotted line) with the 
EPR transmitters used in the single-cell readout (solid lines). The 
lower solid line refers to Q($|l, n), i.e., a single energetic TMSV 
state p epr (l, 1, n), while the upper solid line refers to Q($|oo, n), 
i.e., the optimal EPR transmitter p epr (oo,oo,n) corresponding to 
infinite copies of TMSV states with vanishing signal energy. Lower 
panel. As in the upper panel, except that now we compare C c ($|n) 
with Q($|l, n/2) and Q($|oo, n/2). Despite we assume a stronger 
energy constraint involving the mean total number of photons in both 
signal and reference modes, the single-cell quantum reading is still 
able to outperform the asymptotic multi-cell classical reading. 

Finally, it is interesting to check if the single-cell quantum 
reading represents a superior readout strategy even if we con- 
sider a stronger energy constraint, for instance if we fix the 
mean total number of photons in both the signal and refer- 
ence modes for each copy of the transmitter. Note that this 
approach has been first considered in Ref. D9I1 for individu- 
ating the optimal thermal probes (i.e., the optimal squeezed 
thermal vacua) for detecting the presence of loss in bosonic 
channels. While this stronger energy constraint does not make 
any difference for the classical reading capacity (since the 
optimal classical transmitter involves signal modes only) it 
clearly affects the EPR transmitters where the mean total en- 
ergy of the TMSV states is split exactly in two between sig- 
nal and reference modes. Thus, imposing this stronger energy 
constraint corresponds to compare C c ($|n) with Q($|l,n/2) 
and Q($|oo, n/2). As shown by the lower panel of Fig. QT| 
we see that the single-cell quantum reading is still able to beat 
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the asymptotic multi-cell classical reading. 



VIII. CONCLUSION 
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In this paper we have extended the basic model of quan- 
tum reading to the optimal and asymptotic multi-cell scenario. 
Here the classical memory is modelled as a large block of cells 
where information is stored by encoding a suitable channel 
codeword (channel encoding). This information is then re- 
trieved by probing the whole memory in a parallel fashion and 
detecting the output via an optimal collective measurement 
(channel discrimination). In this general scenario, we define 
the quantum reading capacity of the memory which is a non- 
trivial quantity to compute under the assumption of physical 
constraints for the decoder. In the case of optical memories, 
where data encoding is realized by bosonic channels, the main 
physical constraint is energetic. This leads to define the quan- 
tum reading capacity of an optical memory as the maximum 
number of bits per cell which can be read by irradiating n 
mean photons per cell. 

Despite the general calculation of this capacity is extremely 
difficult, we are still able to provide non-trivial lower bounds 
in the case of optical memories with binary cells. The first 
lower bound, which we have called classical reading capacity, 
represents the maximum number of bits per cell which can 
be read by classical transmitters. This bound has a remark- 
ably simple analytical formula and can be achieved by using 
a single-mode coherent state transmitter. Besides this result, 
we have also computed other bounds by considering particu- 
lar kinds of nonclassical transmitters, including the ones con- 
structed with TMSV states (EPR transmitters), NOON states 
and Fock states. Then we have shown that, in the regime of 
few photons and high reflectivities, these nonclassical trans- 
mitters are able to outperform any classical transmitter, thus 
showing the separation between the classical and the quantum 
reading capacities. It is remarkable that using a single-mode 
or two-mode transmitter per cell is already sufficient to beat 
any classical strategy. Furthermore the classical reading ca- 
pacity can be outperformed even if we restrict the EPR trans- 
mitters to the single-cell readout and we adopt the stronger 
energy constraint where the energy of the reference modes is 
also taken into account. 

In conclusion, our study considers the optimal multi-cell 
encoding for classical memories where we fully extends the 
advantages of quantum reading, i.e., the readout by nonclas- 
sical transmitters. These advantages are particularly evident 
in the regime of few photons with nontrivial consequences for 
the technology of data storage. 
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Appendix A: Miscellaneous proofs 
1. Reduction to pure transmitters 

For the sake of completeness, we show here that the max- 
imization in Eq. (|4oT i can be restricted to pure transmitters. 
This is a trivial consequence of the convexity of the Holevo 
information. 

Let us consider a classical memory with marginal cell <I> 
which is read by an arbitrary transmitter with s signals and r 
references 

p(s,r)= dyp v ipy, ip y = \ip v )(ipy\, (Al) 

where p y ^ and J dy p y = 1 . Then, the conditional Holevo 
information obeys the inequality 



x[$|p(s,r)] < / dy p y x($\il>v) 



(A2) 



In order to prove Eq. dA2l ). let us consider an auxiliary sys- 
tem associated with the variable y. We denote by {\y}} an 
orthonormal basis of this system. Then, we can express the 
transmitter as 



(A3) 



p(s, r) = Ti^ dypyipy® \y)(y\ 



Since the Holevo information cannot increase under partial 
trace, then we have 



(A4) 



X [<£\p(s,r)} ^ X *| / dypyip v ®\y}(y\ 



dyp y x(®\ipy) 



(A5) 



Thus, for any input transmitter p(s, r) we can always choose 
a pure transmitter tf)(s,r) = \ip)(ip\ such that x[ < &|p(s,r-)] sC 
x[$|-0(s, r)). As a result, the maximization in Eq. d46b can be 
restricted to pure transmitters ip(s,r). 



2. Triviality of the unconstrained version of the capacity 

Here we provide a simple sketched proof showing that 
unconstrained quantum reading capacity simply equals the 
whole data stored in the marginal cell of the memory. 

Let us consider a pure transmitter in the tensor-product 
form ip(s) = tp® s at the input of a memory with marginal 
cell $ = {cj) x , p x }. At the output of the cell the arbitrary state 
is given by 



P,{S) = [0xWO]® S , 



(A6) 
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where <fi x (ip) is the single-copy output state. Since the quan- 
tum channels <fi x are different, for any pair <f> x and (f> x i there is 
at least an input (pure) state ip such that 



F[M4>),<i>x>('l>)]=e<l 



(A7) 



For the sake of simplicity, let us assume that this state ip is 
the same for all the channels, i.e., the Eq. (IA7t holds for any 
x 7^ x'. Then, by exploiting the multiplicativity of the fidelity 
under tensor product states, we get 



F[p x (s),p x >(s)] 



(A8) 



for any x 7^ x' . Now, since this quantity goes to zero for 
s — ;> +00 we have that the multi-copy output states p x (s) 



become asymptotically orthogonal. This implies that 

X[*IV(*)] = x({p*(s),P*}) -► H(X) . (A9) 

The proof can be easily extended to the weakest case where 
Eq. (!A7b holds for different input states ipi where i = 
0, • • • , k — 1 for a fc-ary variable X. In this general case, 
for high values of s 3> k, we consider input states 

v>( s ) = vr° ® ■ • • ® ^r 1 

where so H + Sfc-i = s. It is easy to check that the output 

states become asymptotically orthogonal. 
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